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Abstract 

In the chaotic Lorenz system, Chen system and Rossler system, then' equihbria are unstable 
and the number of the equihbria are no more than three. This paper shows how to construct 
some simple chaotic systems that can have any preassigned number of equilibria. First, a chaotic 
system with no equilibrium is presented and discussed. Then, a methodology is presented by 
adding symmetry to a new chaotic system with only one stable equilibrium, to show that chaotic 
systems with any preassigned number of equilibria can be generated. By adjusting the only 
parameter in these systems, one can further control the stability of their equilibria. This result 
reveals an intrinsic relationship of the global dynamical behaviors with the number and stability 
of the equilibria of a chaotic system. 
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Thus far, one typically calculates the number and determines the corresponding 
stability of the equilibria of a given chaotic system, but seldom considers how 
to generate a desirable chaotic system with a preassigned number and preferred 
stability of equilibria. This article shows that, for 3D autonomous systems, 
the latter can be done effectively in a fairly simple way. This helps better 
understand the intrinsic relationship between the equilibria and dynamics of a 
chaotic system. 



I. INTRODUCTION 



In chaos theory, it is important to study the stability of the equilibria of an autonomous 
dynamical system. For a dynamical system described by a set of autonomous ordinary 
differential equations (ODEs), x = /(x), x G -R" , if /(xg) = has real solution then Xg is 
called the equilibrium of this dynamical system. An equilibrium is said to be hyperbolic if all 
eigenvalues of the system Jacobian matrix have non-zero real parts. A hyperbolic equilibrium 
for three-dimensional (3D) autonomous system can be a node, saddle, node-focus or saddle- 
focus. For 3D autonomous hyperbolic type of dynamical systems, a commonly accepted 
criterion for proving the existence of chaos is due to Si'lnikov jl-5|. 

It has also been noticed that although most chaotic systems are of hyperbolic type, there 
are still many others that are not so. For non-hyperbolic type of chaotic systems, they 
usually do not have saddle-focus equilibria, such as those found by Sprott [6|-[9|. 



The well-known Lorenz system 



10|] and also Chen system 
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12 1 both have two unstable 



saddle-foci and one unstable node. They can generate a two-wing butterfly-shaped chaotic 
attractor. Near the center of the two wings, there lies one an unstable saddle-focus. In 
recent years, some chaotic systems are found which can generate three-wing, four-wing, 
and even multi-wing attractors. Observe that, typically, for those symmetrical four- wing 
attractors, near the center of their wings there also lies one unstable saddle-focus. This 
common feature may imply that the number of equilibria basically determines the shape of 
a multi-wing attractor. Therefore, it is interesting to ask: Is it possible to generate a chaotic 
system with an arbitrarily preassigned number of equilibria? Is the number of equilibria 
always determinate the shape of an attractor? 

Furthermore, regarding the stability of the equilibria, recall that recently Yang and Chen 
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found a chaotic system with one saddle and two stable node- foci 13|, and an unusual 3D 



autonomous quadratic Lorenz-like chaotic system with only two stable node-foci 



IJ]. More- 



over, Wang and Chen found an interesting chaotic system with only one stable node-focus 



15| . Thus, another interesting question is whether it is possible for a chaotic system to have 



two, or three, or even an arbitrarily large number of stable/unstable equilibria? 

This paper attempts to answer these questions by showing a novel example of a chaotic 
system with no equilibrium, and several other examples of chaotic systems which have 
any preassigned number of equilibria. Meanwhile it is shown that, by adjusting the only 
parameter in these systems, one can also control the stability of these equilibria. 



II. CHAOTIC SYSTEM WITH NO EQUILIBRIUM 

First, a 3D chaotic autonomous system system is introduced: 

x = y 

< y = z (1) 
z = —y + 3y^ — x"^ — xz — a , 

When a > 0, this system has two symmetrical equilibria: (-y/a, 0,0) and (— 1/0,0, 0), 
as shown by Fig. [H When a = 0, these two symmetrical equilibria merge into one, the 
origin (0,0,0). When a < 0, there is no equilibrium in this system, but still, the system 
can generate a chaotic attractor, as shown in Fig. El The largest Lyapunov exponent with 
respect to the parameter a is shown in Fig. [3l which convincingly implies that the system 
is chaotic. 



III. A MODIFIED SPROTT E SYSTEM WITH ONE STABLE EQUILIBRIUM 

Now, a chaotic system with only one equilibrium, a stable node-focus, is introduced. 



which was reported in 



15|- 



This system was found by modifying the Sprott E system, as follows: 



X = yz + a 
< y = x"^ — y 



(2) 
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FIG. 1: (Color online) The chaotic attractor of system ([T]), which has two symmetrical equilibria 
when a = 0.01: 3D views on the x-y plane, x-z plane and y-z plane. 

When a = 0, it is the Sprott E system [8|; when a 7^ 0, however, the stabihty of the single 
equihbrium is fundamentally different. 

Specifically, when a > 0, system ([2]) possesses only one stable equilibrium: 



P{xE,yE,ZE) 



1 1 

4' 16' 



-16 a 



(3) 



Some numerical calculation results are shown in Table 1. 

Interestingly, this system ([2]) can generate a one-scroll chaotic attractor, as shown in Fig. 

SI 

In the following, this system ([2]) is further modified by imposing some kind of symmetry 
onto it, to have different numbers of equilibria while keeping this system chaotic. 
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FIG. 2: (Color online) The chaotic attractor of system[Tl which has no equilibrium when a = —0.05; 
3D views on the x-y plane, x-z plane and y-z plane. 
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FIG. 3: (Color online) The largest Lyapunov exponent of system ([1]) with respect to the parameter 
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Three-dimensional view 
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FIG. 4: (Color online) The chaotic attractor of system ([2]), when a = 0.006; 3D views on the x-y 
plane, x-z plane and y-z plane. 

TABLE I: Jacobian eigenvalues of system ([2]) 



a 


Jacobian eigenvalues 




0.006 
0.022 


-l,±0.5i 
-0.9607, -0.0197 ±0.5098i 
-0.8458,-0.0771 ± 0.5382i 
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IV. CHAOTIC SYSTEM WITH TWO EQUILIBRIA 



Rewrite system ([2]) in terms of u, v, w as follows: 



u = vw + a 

V = v? — V (4) 
w = 1 — Au . 

From Fig. HJ one can see that the y-axis does not intersect with the attractor in system 
). So, one may try to add a y-axis rotation symmetry to this system, as detailed below. 
Consider the following simple coordinate transformation: 



u = X 
v = y 
w = 



(5) 



2xz . 



This transformation can add a y-axis rotation symmetry, ]Rj,(7r), to the original system, 
because for each {u, v, w) there are two points (±x, ±y, ±z) corresponding to (m, v, w). 
After the above transformation, the system becomes 



• 1 z+2yx'^z+xa-4:x'^z+4:Z^ 

^ ~ 2 J^+P 



y 

z 



[X 



(6) 



1 2yxz^+za~4xz'^ —x+4x^ 
'2 ^iPT^ 



The new system ([6]) possesses two symmetrical equilibria, which are independent of the 
parameter a: -Pl(|, jq, 0) and Pl(— |, Jqi^). 

System is not globally but only locally topologically equivalent to the original system 
([2]), however. Yet one can control the stability of these equilibria by adjusting the parameter 
a, so that the stability remains the same as the original system ([2]) which, when a < are 
unstable and when a > are stable. 

By linearizing system ([2]) at -Pl(|, ]^)0), one obtains the Jacobian 



J 



o 



16 



-2a 
i -1 
-4 -2a 



(7) 



whose characteristic equation is 

det(A/ - J\o) = + (1 + 4a)A2 + (^a'^ + 4a + ^) A + 40^ + ^ = 
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FIG. 5: (Color online) The new two-petal chaotic attractor with stable equilibria when a = 0.003. 
which yields 

Ai = -1 < 0, 
A2 = -2a + 0.5z, 
A3 = -2a-0.5z. 

System ([6]) can generate a symmetrical two-petal chaotic attractor, as shown in Fig. [5] 
and Fig. [6l respectively. 

Numerical calculation of the largest Lyapunov exponent of the the system indicates the 
existence of chaos for some particular values of parameter a, as shown in Fig. [71 
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FIG. 6: (Color online) The new two-petal chaotic attractor with unstable equilibria when a = 
-0.01. 




FIG. 7: (Color online) The largest Lyapunov exponent of system Q with respect to parameter a. 
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V. CHAOTIC SYSTEM WITH THREE EQUILIBRIA 

Similarly, consider the following transformation: 

u = — Sxz"^ 
< v = y (8) 
w = Sx'^z — z^ . 

It can add a y-axis rotation symmetry, ]Rj,(|7r), to the original system, because for each 
{u,v,w) there are three symmetrical points corresponding to it. 
After the above transformation, the system becomes 

• 1 3x'^zy—4x^z^y+x^a—8x^z+2zx+24x'^z''^+z^y~z^a 

^ ~ 3 2x^z^+x^+z-^ 

< y= (x^ - Sxz^f - y (9) 

• _ 1 6 z^x^y-2 z'^xy+2 zxa+4 x^ -x'^ -l& z^x^+z^ +12 z'^x 
^ ~ 3 2i^^Ha;4+24 " 

This system possesses three symmetrical equilibria, which are dependent on the parameter 
a. The analytical expression are too long to write out here, so only some numerical results 
are shown in Table [TTl 



TABLE II: Equilibria and Jacobian eigenvalues of system ([9]) 





a 


Equilibria 


Jacobian eigenvalues 


Unstable case 


-0.01 


Pl=(0.6550,0.0625,0.1258) 

P2=(-0.2186,0.0625,-0.6300) 

P3=(-0.4365,0.0625,0.5044) 


-1.0617,0.0308 ±0.4843i 


Stable case 


0.01 


Pl=(0.6550,0.0625,-0.1258) 
P2=(-0.2186,0.0625,0.6300) 
P3=(-0.4365,0.0625,-0.5044) 


-0.9334, -0.0333 ± 0.5165i 



The system has a symmetrical three-petal chaotic attractor, as shown in Fig. |8]and Fig. 
[9l respectively. 

Numerical calculation of the largest Lyapunov exponent indicates the existence of chaos 
for some particular values of parameter a, as shown in Fig. [101 
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FIG. 8: (Color online) Chaotic attractor of system ([9]) with three stable symmetrical equilibria 
when a = 0.01. 

VI. CHAOTIC SYSTEM WITH ANY NUMBER OF EQUILIBRIA 

Theoretically, one can use the transform (x + iz)^ = {u + iv) to obtain a new chaotic 
system with n equilibria. The following is such an example with five equilibria. 
Consider the following transformation: 

/■ 

u = — 10x^2;^ + bxz"^ 
< v = y (10) 
w = bx^z — IQx'^z^ + z^ . 

It can add a y-axis rotation symmetry, Rj^(|7r), to the original system. The new equations 
are too long to write out, so are omitted here. The symmetrical attractor generated with 
parameter a = is shown in in Fig. [11] 
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FIG. 9: (Color online) Chaotic attractor of system ([9]) with three unstable symmetrical equilibria 
when a = —0.01. 
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FIG. 10: (Color online) The largest Lyapunov exponent of system Q with respect to the parameter 
a. 
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FIG. 11: (Color online) Chaotic attractor of the new system with five symmetrical equilibria when 
a = 0. 

VII. DISCUSSIONS 

The Hartman-Grobman Theorem is an important theorem in ODE systems theory. It is 
about the local behavior of an autonomous dynamical system in the neighborhood of a hy- 
perbolic equilibrium, stating that the behavior of the dynamical system near the hyperbolic 
equilibrium is qualitatively the same as (i.e., topologically equivalent to) the behavior of its 
linearization near this equilibrium point. 

Notice, however, that the new systems discussed in this paper have chaos, which is a 
global behavior, although such a system has only one hyperbolic equilibrium point. In 
other words, all system flows locally converge to the stable equilibrium, but they are chaotic 
globally. This interesting phenomenon is worth further studying, both theoretically and 
experimentally, to further reveal the intrinsic relationship between the local stability of an 
equilibrium and the global complex dynamical behaviors of a chaotic system. 
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